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Method and apparatus for obtaining heartbeat measurements in an EGG waveform 



(54 

(57) A method and apparatus lor obtaining heart- 
beat measurements obtains ECG data from a plurality 
of ECG waveforms, which are in turn obtained from sig- 
nals received from a plurality of ECG electrodes. QRS 
detection logic (71 ) detects heartbeats in the ECG data. 
Classification logic (74) classifies heartbeats into cate- 
gories based on shape and/or timing. Alignment logic 
(75) aligns the heartbeats. Representative heartbeat 
creation logic (77) creates a representative heartbeat 
from the aligned heartbeats. Measurement logic (78) 
measures various aspects of the representative heart- 



beat. This logic analyzes the ECG waveforms to deter- 
mine an earliest QRS onset and latest QRS offset, and 
uses these values to perform a variety of measure- 
ments. This results in robust measurements even in 
very noisy environments. The representative heartbeat 
is displayed, either alone or with heart rate an(dl/or other 
measurement information, to the cardiologist or medical 
professional for diagnosis of the condition of the pa- 
tient's heart, such as a diagnosis of coronary artery dis- 
ease, based on finding a depressed ST segment in the 
representative heartbeat of a patient undergoing a 
stress or exercise test. 
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Description 
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Cross Reference to Related Applications 

This application is related to the following commonly assigned patent applications, filed on even date herewith: 



Title 


Serial Number 


Method And Apparatus For Detecting Heartbeats In An ECG Waveform 




Method And Apparatus For Calculating A Heart Rate In An ECG Waveform 




Method And Apparatus For Classifying Heartbeats In An ECG Waveform 




Method And Apparatus For Creating A Representative Heartbeat From An ECG Waveform 
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Field of the Invention 



2S 



30 



This invention relates to the electronics circuitry field. More particularly, this invention is a method and apparatus 
for obtaining heartbeat measurements from an ECG waveform. 

Baclcground of the Invention 

Designers of medical instrumentation such as cardiographs face many difficult challenges in their jobs. The devices 
they design are expected to deliver high quality information about the electrical activity of a patient's heart to a cardi- 
ologist or other medical professional, so that a correct diagnosis of the condition of the patient's heart can be made. 
Unfortunately, the ECG electrodes connected to a patient usually deliver ECG data to a cardiograph that comprise not 
only information showing the electrical activity of the patient's heart, but also electrical noise. This noise can make up 
most of the ECG data, and can corrupt and totally overwhelm the portion of the ECG data that contain information 
about the electrical activity of a patient's heart. This problem is especially acute in hostile environments, such as a 
patient undergoing a stress or exercise test, where the noise can be quite extreme. Unless the medical instrumentation 
designers are successful at designing a medical instrument that analyzes this ECG data to eliminate or reduce the 
effects of this noise, the cardiologist or other medical professional will find it difficult, if not impossible, to obtain infor- 
mation about a patient's heart, such as measurements of a representative heartbeat, useable in arriving at a correct 
diagnosis of the condition of the patient's heart. 

Summary of the Invention 

A method and apparatus for obtaining heartbeat measurements obtains ECG data from a plurality of ECG wave- 
forms, which are in turn obtained from signals received from a plurality of ECG electrodes. QRS detection logic detects 
heartbeats in the ECG data. Classification logic classifies heartbeats into categories based on shape and/or timing. 
Alignment logic aligns the heartbeats. Representative heartbeat creation logic creates a representative heartbeat from 
the aligned heartbeats. Measurement logic measures various aspects of the representative heartbeat. This logic an- 
alyzes the ECG waveforms to determine an earliest QRS onset and latest QRS offset, and uses these values to perform 
a variety of measurements. This results in robust measurements even in very noisy environments. The representative 
heartbeat is displayed, either alone or with heart rate and/or other measurement information, to the cardiologist or 
medical prof essbnal for diagnosis of the condition of the patient's heart, such as a diagnosis of coronary artery disease, 
based on finding a depressed ST segment in the representative heartbeat of a patient undergoing a stress or exercise 
test. 
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Description of the Drawings 

Fig. 1 shows a block diagram of the medical device of the preferred embodiment of the invention. 
Fig. 2 shows a block diagram of the medical device of Fig. 1 in more detail. 
Fig. 3 shows the processor of the computing unit of the medical device of Fig. 2 in more detail. 
Figs. 4 and 12 show flowcharts of the operation of the cardiograph of the preferred embodiment of the Invention. 
Fig. 5 shows a high level flowchart of the operation of the computing unit of the preferred embodiment of the 
invention. 

Fig. 6 shows a flowchart of the operation of the QRS detection logic of the preferred embodiment of the Invention. 
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Fig. 7 shows a flowchart of the operation of the heart rate calculation logic of the preferred embodiment of the 
invention. 

Fig. 8 shows a flowchart of the operation of the classification logic of the preferred embodiment of the invention. 
Fig. 9 shows a flowchart of the operation of the alignment logic of the preferred embodiment of the invention. 
s Fig. 10 shows a flowchart of the operation of the representative heartbeat creation logic of the preferred embodi- 

ment of the invention. 

Fig. 1 1 shows a flowchart of the operation of the measurements logic of the preferred embodiment of the invention. 
Fig. 1 3 shows a graph of three exemplary ECG waveforms used by the QRS detection logic of the preferred 
embodiment of the invention. 

10 Fig. 14 shows a graph of an exemplary activity function used by the QRS detection logic of the preferred embod- 

iment of the invention. 

Fig. 15 shows a graph of exemplary classified heartbeats. 

Fig. 16 shows a graph of exemplary heartbeats being aligned by the alignment logic of the preferred embodiment 

of the invention. 

7^ Fig. 17 shows a graph of exemplary aligned heartbeats being timesticed by the representative heartbeat creation 

logic of the preferred embodiment of the invention. 

Fig. 18 shows a printout or display of a representative heartbeat without measurements. 
Fig. 1 9 shows a printout or display of a representative heartbeat with measurements. 

20 Detailed Description of the Preferred Embodiment 

I. Overview 

Fig. 1 shows a block diagram of the medical device of the preferred embodiment of the invention. Medical device 
2S 10 comprises acquisition unit 20, electrodes 25, cardiograph 40, and computing unit 60. In the preferred embodiment, 
cardiograph 40 and acquisition unit 20 are separate components of a PageWriter XLi, manufactured by the Hewlett- 
Packard company, modified to execute the flowcharts of Figs. 4 and 12 of the preferred embodiment of the invention. 
Computing unit 60 is a HP Vectra personal computer, suitably programmed to execute the flowcharts of Figs. 5-11 of 
the preferred embodiment of the invention. 
30 Fig. 2 shows a block diagram of medical device 1 0 in more detail. Cardiograph 40 contains acquisition unit interface 

41 , processor 45, printer 47, and computing unit interface 49. Processor 45 executes the flowcharts of Figs. 4 and 12 
of the preferred embodiment of the invention. Computing unit 60 contains cardiograph interface 61, processor 65, 
display 66, input device 67, memory 68, and storage 69. Processor 65 executes the flowcharts of Figs. 5-11 of the 
preferred embodiment of the invention. While Fig. 2 shows medical device 10 as containing discrete components, 
3S those skilled in the art will appreciate that medical device 1 0 could be a single unit that contains each of the components 
shown in Fig. 2, or contain a different number of discrete components, and still fall within the spirit and scope of the 
invention. 

Fig. 3 shows processor 65 of computing unit 60 of medical device 10 in more detail. Processor 65 contains QRS 
detection logic 71, heart rate calculation logic 73, classification logic 74, alignment logic 75, representative heartbeat 
40 creation logic 77, and measurement logic 78. In the preferred embodiment, each of these logic blocks is performed by 
software written to perform the functions of relevant portions of the flowcharts shown in Figs. 5-11, and this software 
is executed by processor 65. Alternatively, some or all of logic blocks 71-78 could be special purpose hardware, such 
as contained in an application specific integrated circuit, designed to perform functions of relevant portions of the 
flowcharts shown in Figs. 5-11. 

"^s Fig. 4 shows a flowchart of the operation of cardiograph 40 of the preferred embodiment of the invention. In block 

101 . ECG signals are received'from electrodes 25 of acquisition unit 20. In the preferred embodiment, these signals 
are digital signals sampled at a high sampling rate. Block 103 decimates and fitters the sampled ECG signals outside 
of a predetermined bandwidth. In the preferred embodiment, the predetermined bandwidth is .01 Hz to 150 Hz, and 
the decimation process reduces the number of samples to one eighth of the number of original samples. Block 105 

50 makes twelve ECG leads from the ten electrodes in a conventional manner. The signals contained on the twelve ECG 
leads will be referred to herein. as "ECG waveforms", and the information contained thereon will be referred to herein 
as 'ECG data". Those skilled in the art will appreciate that the actual number of electrodes or leads may be different 
than discussed above and still fall within the spirit and scope of the invention. Block 108 sends the ECG data on the 
ECG waveforms to computing unit 60. The flowchart ends in block 109. 

^ Fig. 5 shows a high level flowchart of the operation of computing unit 60 of the preferred embodiment of the 

invention. Block 201 receives the ECG waveforms containing the ECG data from cardiograph 40. Block 203 fon^/ard 
fitters the ECG data. In the preferred embodiment, this filter is a high pass filter used as part of a forward/reverse 
filtering scheme to remove baseline wander while preserving low frequency information in the ECG data. Block 300 
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calls a subroutine that detects the heartbeats (i.e. QRS complexes) in the ECG waveforms. This subroutine calculates 
an activity function from a subset of the ECG waveforms determined to be least noisy, and uses this activity function 
to search for heartbeats. This allows for true heartbeats to be detected while discarding false "noise" beats. The op- 
eration of this subroutine will be described in more detail later in conjunction with the discussion of Fig. 6. 

s Block 400 calls a subroutine that calculates the patient's heart rate. This logic determines the intervals between 

the heartbeats, discards a percentage of the shortest and longest intervals, and averages the remaining intervals to 
arrive at the patient's heart rate. This results in a robust calculation of the heart rate even in the presence of noise 
falsely detected as heartbeats and missed beats common in noisy environments. The operation of this subroutine will 
be described in more detail later in conjunction with the discussion of Fig. 7. 

10 Block 500 calls a subroutine that classifies heartbeats. This classification is done by comparing each heartbeat 

against a group of templates corresponding to one or more heartbeat classifications. The templates are updated to 
track changes in the morphology of the heartbeats. The operation of this subroutine will be described in more detail 
later in conjunction with the discussion of Fig. 8. 

Block 205 reverse filters the ECG data. In the preferred embodiment, this filter is a high pass filter used as part of 

IS a forward/reverse filtering scheme to remove baseline wander while preserving low frequency Information in the ECG 
data. 

Block 600 calls a subroutine that aligns heartbeats prior to representative heartbeat creation. This logic slides the 
heartbeats across an alignment template heartbeat to calculate when the heartbeats are aligned, and performs ad- 
justments to reduce the effects of noise or jitter on the different ECG waveforms. The operation of this subroutine will 

20 be described in more detail later in conjunction with the discussion of Fig. 9. 

Block 700 calls a subroutine that creates a representative heartbeat from the aligned heartbeats. This logic time 
slices through the aligned heartbeats, discarding a percentage of the smallest and largest magnitudes of the aligned 
heartbeats at each instance of time and averaging the remaining magnitudes to produce a representative heartbeat. 
This trimmed averaging technique results in a high quality representative beat, since samples from noise and misclas- 

2S sified beats are discarded. The operation of this subroutine will be described in more detail later in conjunction with 
the discussion of Fig. 10. 

Block 800 calls a subroutine that measures various aspects of a representative heartbeat. This logic analyzes the 
representative heartbeats from a group of ECG waveforms to determine an earliest QRS onset and latest QRS offset, 
and uses these values to perform a variety of measurements. This results In robust measurements even In very noisy 
30 environments. The operation of this subroutine will be described In more detail later in conjunction with the discussion 
of Fig. 11. 

Block 210 displays the representative heartbeat created by subroutine 700 and, optionally, the measurements 
obtained by subroutines 800 and 400, on display 66 of computing unit 60 (Fig. 2). Examples of these displays are 
shown in Figs. 18 and 1 9. Block 220 sends the representative heartbeat and measurements for each ECG waveform, 

3S including the heart rate measurement calculated by subroutine 400, back to cardiograph 40. Cardiograph 40 processes 
this information in accordance with the flowchart of Fig. 1 2. The flowchart ends in block 249. 

Fig. 1 2 shows how cardiograph 40 processes the information received from computing unit 60. Block 1 50 receives 
the representative heartbeat and measurements, including the heart rate measurement, sent by block 220 of Fig. 5. 
Block 190 prints the representative heartbeat created by subroutine 700 and, optionally the measurements obtained 

40 by subroutines 800 and 400, on printer 47 of cardiograph 40 (Fig. 2). Examples of these printouts are shown in Figs. 
18 and 19. 

II. QRS Detection 

"^s Fig. 6 shows a flowchart of the operation of subroutine 300, performed by QRS detection logic 71 of computing 

unit 60 of the preferred embodiment of the invention. Block 302 performs an initialization process that learns about the 
ECG data. More specifically, the first time through the subroutine, this process analyzes the first few seconds of ECG 
data to determine a preliminary interval between peaks and the magnitude of an average peak. During routine process- 
ing (i.e. subsequent times through the subroutine) block 302 continues to update the information obtained the first time 

50 through the subroutine. Block 301 receives ECG data from three ECG waveforms out of the twelve ECG wavefonns 
received by computing unit 60 In bbck 201 (Fig. 4). A graph of three exemplary ECG waveforms is shown in Fig. 13. 
In the preferred embodiment, block 301 selects these three ECG waveforms as the three ECG waveforms that have 
optimal noise characteristics. This is done by continuously calculating the signal quality on the twelve ECG waveforms 
and ranking these waveforms from highest to lowest signal quality. 

ss Block 303 calculates an activity function out of the three leads of ECG data. An activity function Is a signal math- 

ematically derived from the ECG data which emphasizes characteristics of the heartbeat while minimizing the influence 
of noise to enable more accurate heartbeat detection. In the preferred embodiment, the activity function is created by 
decimating the ECG data by one half, filtering the data using a bandpass filter and taking the absolute value of the first 
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difference. The absolute first differences fronn the three ECG waveforms are then summed together, thresholded and 
smoothed to improved noise performance. Noise statistics (i.e., signal quality ) across all twelve ECG waveforms are 
calculated and updated. A graph of an exemplary activity function is shown In Fig. 14. 

Block 305 checks the noise statistics to see If unacceptably high noise was detected. If so, block 306 issues a 
s "high noise" message, which is displayed on display 66 of computing unit 60 (Fig. 2). This error causes the subroutine 
to terminate abnormally by ending in block 398. In the absence of high noise, block 305 is answered negatively Block 
309 looks to see if It Is to detect another heartbeat. If not, the subroutine returns In block 399 to block 400 of Fig. 5. 

If block 309 is answered affirmatively, subroutine 300 performs up to three different types of searches to find each 
heartbeat. The first search is an on-time search. This search is nornnally used to detect heartbeats that appear within 
10 a small window of their expected time. If the on-time search fails to detect a heartbeat, a modified on-time search is 
used. The modified on-time search is similar in operation to the on-time search, but can detect heartbeats in low noise 
environments that the on-time search might miss, such as early beats. If both the on-time search and the modified on- 
time search fail to detect a heartbeat, a late beat search Is performed. This search detects heartbeats that appear later 
than their expected time. 

IS The operation of the on-time search will now be discussed In more detail. Block 310 performs an on-time search 

of the activity function for QRS complexes in the EGG data. In the preferred embodiment, block 310 searches over 
1 1 5% of the current interval length between heartbeats (as determined in block 302) with a threshold that starts at 80% 
of the average peak value of the activity function. The threshold is then linearly decreased overtime to 40% at the end 
of the search window. Any local maxima discovered above this linearly decreasing threshold is considered a "peak". 

20 Block 315 checks to see If any peaks were detected. If so, block 320 selects the true heartbeat (i.e., QRS peak) 

from the detected peaks by looking at the times each peak occurred. Those skilled in the art will appreciate that the 
above-described search could detect multiple peaks in a noisy environment, all but one of whk^h would be noise. The 
peak that occurred closest to the time the next heartbeat was expected is considered a true heartbeat, and timing 
information about this heartbeat is stored in storage 69 of computing unit 60 (Fig. 2). Flow of control loops back to 

2S block 309 to look for another beat to detect. 

Referring again to Fig. 6, the operation of the modified on-time search will now be discussed. If block 315 is 
answered negatively, block 330 checks to see if the noise level Is low. If so, block 335 searches the activity function 
by using a modified on-time search. In the preferred embodiment, this search is performed by using a constant lower 
threshold, such as 54% of the average peak activity function value, over 115% of the current inten^al length. Block 335 

30 checks to see if any peaks were detected. If so, block 338 selects the true heartbeat from the detected peaks by 
selecting the peak that occurred closest to the time the next heartbeat was expected. Infornnation about the heartbeat 
selected in block 338, such as beat timing Information, is stored in storage 69 of computing unit 60 (Fig. 2). Flow of 
control loops back to block 309 to look for another beat to detect. 

The operation of the late beat search will now be discussed. If either blocks 330 or 335 are answered negatively 

3S block 350 performs a late beat search. In the preferred embodiment, this block searches the activity function over three 
R-R interval lengths using a linearly decreasing threshold, as was done with the on-time search. Block 355 checks to 
see if any peaks were detected. If so, block 358 selects the true heartbeat from the detected peaks by selecting the 
first peak that it finds, information about the heartbeat selected in block 358, such as beat timing information, is stored 
in storage 69 of computing unit 60 (Fig. 2). Flow of control loops back to block 309 to look for another beat to detect. 

"^0 If block 355 Is answered negatively block 370 Issues a "detector failed" message that is displayed on display 66 of 
computing unit 60. Since no beats were detected, this causes an abnormal termination of the subroutine in block 396. 

III. Heart Rate Calculation 

Fig. 7 shows a flowchart of the operation of subroutine 400, performed by the heart rate calculation logic 73 of 
computing unit 60 of the preferred embodiment of the invention. Block 401 sets a counter called RR_ctr = 0. Block 403 
reads and stores the R-R interval between the first heartbeat detected in the EGG waveform and the second heartbeat 
detected in the EGG waveform. In the preferred embodiment, this detection is done using information stored by QRS 
detection logic 71 In subroutine 300, although conventional methods of detecting R-R intervals In an ECG waveform 
50 could also be used. Block 405 increments RR_ctr. Block 410 checks to see If the counter is less than a maximum 
counter value and there are additbnal heartbeats still available for analysis in the EGG waveform. If both of these 
conditions are true, flow of control loops back to block 403 until one condition is no longer true. An alternate embodiment 
has been contemplated where block 403 is a timer and where block 410 checks to see if a maximum time has elapsed. 
For example, if a maximum time was set for 10 seconds, only heartbeats occurring In the most recent 10 second period 
S5 of time are used to calculate the heart rate. 

When block 410 Is eventually answered negatively block 420 checks to make sure at least a minimum number of 
heartbeats were analyzed by the loop made up of blocks 403-410. 

If block 420 is answered affirmatively, block 425 sorts the R-R intervals from shortest to longest. Block 430 then 
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discards a percentage of the shortest and longest R-R intervals. In a noisy environment, the QRS detector may incor- 
rectly detect noise as a heartbeat, and may incorrectly miss a real heartbeat. These errors result in incorrect R-R 
inten/als, both too short and too long. The trimmed average done here results in a robust and accurate heart rate 
calculation even In the presence of false detections and missed beats. With low noise and/or arrhythmias, this trimmed 

5 average technique also produces an accurate calculation of the heart rate. 

In the preferred embodiment, block 430 discards 25% of the shortest and 25% of the longest R-R intervals . al- 
though different values could be used. Block 435 then averages the remaining R-R intervals. Block 440 converts this 
average R-R interval to a heart rate. Block 445 smoothes the heart rate determined by block 440 by averaging it with 
a predetermined number of past heart rates. In the preferred embodiment, block 445 averages the current heart rate 

10 with the past two heart rates. In' any event, the heart rate determined by block 445 (or by block 440, if the smoothing 
step of block 445 is not desired) Is stored by block 450 in storage 69 of computing unit 60 (Fig. 2). The subroutine 
returns in block 499 to block 500 of Fig. 5. 

If block 420 is answered negatively, block 460 simply computes the average R-R interval of the small number of 
R-R intervals that were read and stored in block 403. This average R-R interval is converted to a heart rate in block 

IS 440, and the heart rate is stored by block 450 in storage 69 of computing unit 60 (Fig. 2). As before, the subroutine 
retums in block 499 to block 500 of Fig. 5. 

IV. Heartbeat Classification 

20 Fig. 8 shows a flowchart of the operation subroutine 500, performed by classification logic 74 of computing unit 

60 of the preferred embodiment of the invention. In the preferred embodiment, subroutine 500 is used to classify beats 
as "D" (dominant), "V (ventricular ectopic), "S" (supraventricular ectopic), or "Q" (questionable), although other clas- 
sifications could be used. 

Block 501 checks to see if there are more heartbeats to classify. If so, Block 502 gets beat timing information for 
25 the next heartbeat. In the preferred embodiment, this information is obtained from QRS detection logic 71 in a manner 
that has already been discussed, although conventional methods of obtaining this information could also be used. 
Block 503 normalizes the portion of the activity function (from QRS detection logic 71 or conventional means) around 
the beat to classify. Block 505 performs timing and physiologic measurements on both the activity function and on the 
beat to classify. These measurements are used in blocks 510 and 550 to assist in beat classification, as will be discussed 
30 later. 

Block 510 checks to see if all measurements made by block 505 are within physiologic limits (e.g. within a prede- 
termined width and height). If not, the beat is classified as questionable ("Q"), due to noise, and flow of control loops 
back to block 501 to see if there are more beats to classify. If block 510 is answered affirmatively, block 515 sets a 
template counter to 1 , initializing a series of steps that compares the beat to classify with one or more templates. Block 

3S 518 checks to see if a counter called Numjemplates = 0. If so, no templates have yet been created for this ECG 
waveform, and flow of control skips down to blocks 539 and 540 to create a new template by saving this beat as the 
first template. Block 542 then uses timing and physiologic infornnation to classify this beat and template as either 'D' 
(Dominant), "V (ventricular ectopic), "S" (supraventricular ectopic), or "Q" (questionable). Most commonly, this beat 
will be classified as D and the template classified as D/S, for "Dominant/supraventricular ectopic", since the vast majority 

40 of beats classified will be classified this way, and since both D and S beats have the same morphology and thus would 
match the same template, but vary by timing information - the S beat being earlier than the D beat. The method of 
distinguishing between D and S beats is performed by blocks 545 to 558, as will be discussed in more detail later 
Those skilled in the art will appreciate that more than one template can be classified as "D/S", since dominant heartbeats 
can have more than one unique morphology. Block 542 increments Num_Template to one to indicate one stored tem- 

45 plate. Flow of control moves to block 545, the operation of which will be discussed later. 

Once at least one template Is created, block 518 is answered negatively, and block 520 slides the beat to classify 
across the first, stationary template. As previously discussed, the first template usually corresponds to a first heartbeat 
classification of D/S, meaning the template for dominant and supraventricular waveforms. As block 520 slides this beat 
across the first template, the minimum area difference between the beat to classify and the first template is calculated. 

50 Block 525 asks if this minimum area difference is less than a threshold. If it is, the beat matches the first (D/S) template, 
and block 530 averages the new beat data into the template it matched. In the preferred embodiment, a weighted 
average is used, where the existing template is given more weight in the averaging process than the new beat. Block 
530 also keeps track of the number of times a beat matched this template, as well as when a beat most recently 
matched this template. 

55 If block 525 is answered negatively, block 535 increments the template counter. Block 538 verifies that the maximum 

number of templates to check has not been exceeded, which would indicate all templates have been checked. If block 
538 is answered affirmatively, flow of control loops back to blocks 518 and 520, where the beat is slid across the second 
template. The second template exemplarily corresponds to a classification of ventricular ectopic ("V"). Block 520 again 
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computes the minimum area difference, and block 525 again asks if this minimum area difference is less than a thresh- 
old. If it is, the beat matches the second (V) template, and block 530 averages the new beat data into the template it 
matched. If no match was found, the subroutine loops through blocks 535, 538, 518, 520, and 525 until a match is 
found or until block 538 is answered negatively, indicating that all existing templates have been checked for matches. 

s Once a match is found and block 530 averages the new beat data into the template it matched, block 537 verifies 

the classification of the template. As new beats are averaged into the existing templates, it is possible for the classifi- 
cation of the template to change. For example, a template originally classified as "V" may be reclassified to "D" as 
more beats are averaged into It. Block 545 checks to see if the beat matched a D/S template. If so, an additional inquiry 
about the beat must be made before the beat can be classified. This is done in block 550. which asks if the beat was 

10 early. If so, the beat Is classified as supraventricular ("S") in block 555, and flow of control loops back to block 501 to 
look to see if there are more beats to classify. If not, the beat is classified as dominant ("D") in block 558, and flow of 
control loops back to block 501 to look to see if there are more beats to classify. If block 545 determined that the beat 
matched a template other than the D/S template, block 560 classifies the beat as the classification corresponding to 
the template it matched. For example, if the beat matched the ventricular ectopic ("V") template, the beat would be 

IS classified as ventricular ectopic. Flow of control loops back to block 501 , as previously discussed. 

Referring again to block 538, if block 538 is answered negatively, all existing templates have been checked and 
none of them match. Block 539 then asks if Numjemplate is less than Maxjemplate -- a counter indicating the max- 
imum number of templates. If not, block 540 creates a new template for this beat. Block 542 classifies the new template 
and beat and increments Num_templates, as discussed previously. 

20 If block 539 is answered negatively, the maximum number of templates has been reached. Block 541 overwrites 

the template with the least recent update. Like block 542, block 543 classifies the new template, but does not increment 
Num .templates, since the number of templates did not change. Those skilled in the art will appreciate that the number 
of templates actually created can vary, depending on the amount of noise in the environment and whether any ectopic 
beats are detected. 

2S When block 501 determines that there are no more beats to classify, block 590 displays the classified heartbeats 

on display 66 of computing unit 60 (Fig. 2). One such exemplary display is shown in Fig. 15. The subroutine returns 
in block 599 to block 205 of Fig. 5. 



V. Beat Alignment 

30 

Fig. 9 shows a flowchart of the operation of subroutine 600, performed by alignment logic 75 of computing unit 60 
of the preferred embodiment of the invention. Block 601 determines which 3 ECG waveforms are the least noisy. In 
the preferred embodiment, this is done by using the information obtained in block 301 of the Detect QRS Complexes 
in ECG Waveform subroutine 300, although this determination could be made directly in this subroutine by continuously 

36 calculating the signal to noise ratio or other indication of signal quality on the twelve ECG waveforms and ranking these 
waveforms from highest to lowest , or by using some other technique. Block 603 sets a lead counter to look at the first 
of the three least noisy ECG waveforms. Block 605 receives ECG data from the ECG waveform determined by the 
lead counter. After the ECG data is low pass filtered in block 608, block 61 0 determines the dominant template amongst 
the beats to be used to construct the representative heartbeat. In the preferred embodiment, this is done by looking 

40 at the number of times the templates used by the classification logic were matched by a beat, as determined in block 
530 in Fig. 8. By definition, this will be a D/S template. 

Block 620 creates and stores an alignment template in storage 69 of computing unit 60. In the preferred embodi- 
ment, the alignment template is created with a beat which matched the dominant template determined in block 610, 
where a portion of the dominant template around the QRS is normalized . An alternate embodiment has been contem- 
plated where step 610 is skipped, and block 620 creates an alignment template by finding the first beat classified as 
D (via beat classification logic 74 or a conventional method of beat classification) and using this beat as the alignment 
template. 

Block 630 sets a beat counter to 1 . Block 640 gets the next dominant beat for this ECG waveform and normalizes 
a portion of this beat around the QRS complex. For the purposes of this application, beats classified as "D" are referred 

so to herein as "dominant", while beats classified as "V", "S", or "Q" are referred to as "non-dominant". In the preferred 
embodiment, non-dominant beats are excluded from alignment and from the determination of a representative heart- 
beat, since these beats can contaminate the representative heartbeat. 

Block 650 slides this beat across the stationary alignment template, while computing the value of the sum of the 
absolute values of the difference between the beat and the alignment template. This value is referred to herein as the 

ss area difference. The position where the area difference is the minimum (Minimum Area Difference) is the position where 
the beat is best aligned with the alignment template, and this position is saved for this beat in storage 69 of computing 
unit 60. Block 655 increments the beat counter. Block 660 checks to see if the beat counter is greater than or equal to 
the number of beats to align for this ECG waveform. If not. flow of control loops back to block 640 to get the next 
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dominant beat. If so, block 670 incrennents the lead counter 

Block 675 checks to see it the lead counter is greater than or equal to 3 - the nunnber of least noisy ECG waveforms 
selected in block 601 . If block 675 is answered negatively, flow of control loops back to block 605, where the alignment 
process is repeated for the next ECG waveform. If block 675 is answered affirmatively, block 680 corrects and stores 
s each beat time on each ECG waveform with the median value of the alignment times just determined independently 
on the three least noisy leads. This is done to minimize the effect of noise which results in beats on different ECG 
waveforms best aligning at slightly different times (i.e. jitter). The subroutine returns in block 699 to block 700 in Fig. 
5. Fig. 16 shows a new beat being slid across a stationary template. 

10 VI, Representative Beat Creation 

Fig. 10 shows a flowchart of the operation of subroutine 700, performed by representative heartbeat creation logic 
77 of computing unit 60 of the preferred embodiment of the invention. Block 701 determines which beats to use in 
creating the representative heartbeat. In the preferred embodiment, only "dominant" beats, aligned via the beat align* 

IS ment steps discussed above, are used. If more "dominant" beats are available than the number needed to construct 
a representative heartbeat, those with the most similar morphology are used. For example, if two or more dominant 
templates were created in classification subroutine 500, only the beats that match the dominant template that contains 
the most beats will preferably be used.. Block 705 sets a counter which keeps track of which ECG waveform the 
representative heartbeat is being created for to one. 

20 Block 710 checks to see if the lead number counter is greater than or equal to the maximum number of ECG 

waveforms. If not, block 715 reads the aligned "dominant" beats for this ECG waveform. Block 720 sets a time pointer 
to zero. 

Block 725 gets a time slice of data for each aligned beat at the instance of time identified by the time pointer. In 
the preferred embodiment, this data is the magnitude of each of the aligned dominant heartbeats at this moment of 

25 time. Block 730 sorts the magnitudes'for this slice of time from smallest to largest. Block 735 discards a percentage 
of the smallest and largest magnitudes. In a noisy environment, the beat classification subroutine may incorrectly 
classify beats as being dominant. These errors result in misclassified beats being incorrectly included in the aligned 
beats. The trimmed average done here results in a robust and accurate representative heartbeat even in the presence 
of misclassified beats and high noise present on dominant beats. In the preferred embodiment, 33% of the smallest 

30 magnitudes and 33% of the largest magnitudes are discarded, although other values could be used. Block 740 averages 
the remaining magnitudes for this time slice. Block 750 stores the average magnitude for this time slice in a represent- 
ative beat array in storage 69 of computing unit 60. Block 755 increments the time pointer to the next slice of time, and 
block 760 checks to see if the time pointer has reached its maximum value. If not, flow of control loops back to block 
725 to determine the average magnitude of the other slices of time to complete the representative heartbeat array. Fig. 

35 17 shows exemplary aligned beats being time sliced using the process described above. 

When block 760 is answered affirmatively, block 765 fonward and reverse filters the representative heartbeat stored 
in the representative heartbeat array, and stores the result back into the array in block 770. An alternate embodiment 
has been contemplated where this step is skipped. Block 775 increments the ECG waveform counter, and flow of 
control loops back to block 710 to create a representative heartbeat for each of the other ECG waveforms. Once block 

40 710 determines that representative heartbeats have been created and stored for each of the ECG waveforms, block 
780 displays the representative heartbeats on display 66 of computing unit 60. An exemplary display of a representative 
heartbeat is shown In Fig. 18. The subroutine returns in block 799 to block 800 of Fig. 5. 

VII. Measurements 

45 

Fig. 11 shows a flowchart o^ the operation of subroutine 800, performed by measurements logic 78 of computing 
unit 60 of the preferred embodiment of the invention. Block 801 obtains the representative heartbeats for each of the 
ECG waveforms. In the preferred embodiment, this is done by reading the representative heartbeat array stored in 
step 770 of Fig. 10. Alternatively, representative heartbeats created using a different method, including conventionally 

50 known methods for creating representative heartbeats, could be used. 

Block 803 measures the earliest QRS onset and the latest QRS offset. These values are from an activity function 
derived from a subset of the representative heartbeats obtained in block 801. These values are used- for many of the 
measurements that will be made for these representative heartbeats, as will soon be discussed. Block 805 sets a 
counter that keeps track of the ECG waveform for which the representative heartbeat is being measured. Block 810 

55 gets the representative heartbeat for this ECG waveform. Block 81 5 detenmines the isoelectric level of the represent- 
ative heartbeat. In the preferred embodiment, this is the average level of the 16 msec of data prior to the earliest QRS 
onset. Block 820 determines the R wave amplitude for this representative heartbeat. In the preferred embodiment, this 
is the maximum positive value between the earliest QRS onset and the latest QRS offset, with an adjustment made to 
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correct tor elevated ST segments at the latest QRS offset, if the "T".wave is so large that it impacts the determination 
of the R wave amplitude. 

Blocl^ 825 determines the ST level. In the preferred embodiment, this is the average of 10 msec around the user- 
determined ST measurement point of the representative heartbeat. Block 830 determines the ST slope. In the preferred 
s embodiment, this is determined by using a best line fit between the latest QRS offset and the ST measurement point 
of the representative heartbeat. Block 835 determines the ST integral. In the preferred embodiment, this is determined 
by computing the sum of the negative area between the latest QRS offset and the ST measurement point of the rep- 
resentative heartbeat. 

Block 850 updates a measurement confidence flag for each measurement taken. In the preferred embodiment, 
10 historical information and physiologic limits are used to set these measurement flags to either a "low" or "high" confi- 
dence. These confidence flags can be displayed to a user in a variety of ways, including the term "low" or "high" 
displayed next to a measurement, changing the color of the measurement on the display (e.g., green means high, red . 
means low), etc. A "tow" confidence flag would indicate to the cardiologist or other medical professional that a meas- 
urement is not physiologic or has changed in a non-physiologic manner and should be manually reviewed for correct- 
is ness. Block 855 increments the ECG waveform counter. Block 860 checks to see if the ECG waveform counter exceeds 
the maximum number of ECG waveforms. If not, flow of control loops back to block 810 to repeat the process for the 
other ECG waveforms. If so, block 880 displays the measurements on display 66 of computing unit 60 (Fig. 2). An 
exemplary display of these measurements, displayed along with the representative heartbeat, is shown in Fig. 19. A 
cardiologist looking at the representative heartbeat and the measurements shown in Fig. 19 would see that there is a 
20 depressed ST segment, indicating that the patient undergoing a stress test has coronary artery disease. The subroutine 
returns In block 899 to block 21 0 of Fig. 5. 

I 

Claims 

1. A medical device (10) for obtaining heartbeat measurements from an ECG waveform, comprising: 

detection logic (71) for detecting a plurality of heartbeats in said ECG waveform; 

30 classification logic (74) for classifying said plurality of heartbeats as dominant heartbeats and non-dominant 

heartbeats; 

alignment logic (75) for aligning said dominant heartbeats to create aligned dominant heartbeats; 

■ 

3S representative heartbeat creation logic (73) for creating a representative heartbeat from said plurality of aligned 

dominant heartbeats; 

measurement logic (78) for obtaining measurements for said representative heartbeat, said measurement 
logic further comprising: 

40 

means for obtaining a plurality of representative heartbeats corresponding to a plurality of ECG waveforms; 

means for measuring an earliest QRS onset, responsive to said obtaining means; 

45 means for measuring a latest QRS offset, responsive to said obtaining means; and 

for said representative heartbeat, means for determining a maximum R wave amplitude measurement 
between said earliest QRS onset and said latest QRS offset. 

so 2. The medical device of claim 1 , further comprising: 

for said representative heartbeat, means for determining an average ST level measurement. 

3. The medical device of claim 2, further comprising: 

for said representative heartbeat, means for determining an ST slope measurement by using a best line fit between 
^ the latest QRS offset and an ST measurement point of said representative heartbeat. 

4. The medical device of claim 3, further comprising: 

for said representative heartbeat, means for determining an ST Integral measurement by computing the sum of 
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the negative area between said latest QRS offset and said ST measurement point of said representative heartbeat. 

5. The medical device of claim 1 , further comprising: 

means for determining a level of confidence of said maximum R wave amplitude measurement; and 
means for displaying an indication of said level of confidence. 

6. The medical device of claim 2, further comprising: 

means for determining a level of confidence of said ST level measurement; and 
means for displaying an indication of said level of confidence. 

7. The medical device of claim 3, further comprising: 

means for determining a level of confidence of said ST slope measurement; and 
means for displaying an indication of said level of confidence. 

8. The medical device of claim 4, further comprising: 

means for determining a level of confidence of said ST integral measurement; and 
means for displaying an indication of said level of confidence. 

9. The medical device of claim 1 , further comprising: 

for each representative heartbeat, means for determining an isoelectric level by computing an average level for a 
predetermined amount of time prior to said earliest QRS onset. 

10. The medical device of claim 9, further comprising: 

means for determining a level of confidence of said isoelectric level; and 
means for displaying said level of confidence. 

determining a level of confidence of said isoelectric level measurement; and displaying an indication of said 
level of confidence. 



< 



10 




11 




12 



I 

EP 0 776 630 A1 



CO 

1^ 



0) 

"cS 

a: 



b) 
o 



o 



(0 
0) 



JO 
(0 

O 



c 
o 



o 



Q 
CO 

a: 
o 



o 
o 



o 

b) 
o 



0) 

E 
c 

O) 



J 



c 




o 




CO 


o 


o 






'(75 


o 


(/) 




CO 




O 





c 
E 

0) 

0) 

(0 
0) 




00 



b) 
o 



live 


.ogic 


(0 




•eseni 


tBeal 


Repr 


Hearl 



CD 



O 
(/) 
(/) 
CD 
O 
O 

a. 



CO 
d) 

il 



13 



EP 0 776 630 A1 

I 



101 



( start ) 

Receive ECG Signals 
For Acquisition Unit 



103 



Decimate And Filter 
Samples Outside 
Predetermined Bandwidth 



105. 



108 



Make 12 Leads For 
1 0 Electrodes 



I 



Send ECG Data To 
Computing Unit 



Fig. 4 



C End > 



109 



14 



EP 0 776 630 A1 



201 



C Start ) 



Receive 1 2 Lead ECG 
Data From Cardiograph 



203 



300 



Forward Filter Data 

I 



Detect QRS Complexes 
In ECG Waveform (Fig. 6) 



400 



Calculate Heart Rate 
(Fig. 7) 



500 



f 



Classify Beats In ECG 
Waveform (Fig. 8) 



205 




600 



700 



Align Beats In ECG 
Waveform (Fig. 9) 



Create Representative 

Beat From ECG 
Waveform (Fig. 10) 



800 
J- 



Obtain Measurements 
From Representative 
Beat For ECG Waveform 
(Fig. 11) 



220 



Display Representative 
Beat And Measurements 



I 



210 



Send Representative Beat And 
Measurements Back To Cardiograph 



Fig. 5 



C End y 



249 



15 



EP 0 776 630 A1 



300 



Detect QRS Complexes 
In ECG Waveform 



I 



301 J 



Perform Initialization ["302 




Receive 3 Leads Of Data 



T 



398 

Ylnd) 



Make Activity Function | 



303 



issue 
"High Noise" 
Message 



306 




399 
Return) 



16 



EP 0 776 630 A1 



320 
-I 



Select True Beat 
From Detected 

Peaks And Store 
Beat Time 





Perform Modified 
On-Time Search 



338 

N 



Select True Beat 
From Detected Peal<s 
And Store Beat Time 



Any 
Peaks 
Detectei 
? 



335 



No 



350 



f 



Perform Late-Beat Search 



358 

\ 



Select True 
Beat From 
Detected Peaks 
And Store 
Beat-Time 




No 



370 



Issue 
"Detector-Failed" 



396 



Message 
H End ) 



Fig. 6B 



17 



EP 0 776 630 A1 



400 



Calculate Heart Rate 



Set RR CTR = 0 



401 




Fig. 7A 



Compute 
Average 
R-To-R 
Interval 



18 



EP 0 776 630 A1 

I 



425 
> — 



430 



Sort R-To-R Intervals 



Discard Predetermined % of 
Shortest and Longest R-To-R Intervals 



445 



I 



Average Remaining % of 
R-To-R Intervals 



435 



440 



Convert Average R-To-R 
Interval To Heart Rate 



Smooth By Averaging With 
Past Heart Rate(s) 




X Return) 



Fig. 7B 



19 



EP 0 776 630 A1 



500 



Classify Beats in 
ECG Waveform 



502 




590 



Display 
Classified 
Heartbeats 



599 > ^ V 
^Return) 



Get Beat Timing Information 



I 



J503 



Normalize Portion Of QRS 
Detector Activity Function Around 
The Beat To Classify 



1 



Make Measurements On The 
Activity Function And On The 
Beat To Classify 



505 



Are All 



510 



Measurements 
i/Vithin Physiologi 
Limits ? 



Fig. 8 A 



Yes 



20 



EP 0 776 630 A1 



615 




1 SetCNT = T1 




Yes 




Slide The Unclassified Beat Across 
Template [CNT] While Computing 
The Minimum Area Difference 



520 



535 



Increment 
CNT 



No 



. 538 

"^^^As CNT 
Num_Template§; 
? 

No 



^525 

Is The 
Minimum Area 
Difference < 
Threshold 

7 



Yes 



530 
^ 



Average New 
Beat Data 

Into Template 
It Matched 



537 



I 



Verify 
Template 
Classification 



Fig. 8B 




21 



EP 0 776 630 A1 




(D 




Overwrite 
Template 
Updated 
Longest 
Time Ago 


541 




Classify New 
Template 
And Beat 


543 









Classify New 
Template And 
Beat And 
Increment 
Num_Templates 




560 
-J. 



Classify As 
Classification 
Corresponding 
To Template 
Beat Matched 



558 



No 


Classify 






AsD 





Fig. 8C 



22 



EP 0 776 630 A1 



60QJ 



Align Beats In ECG Waveform 



I 



601s 



Determine Identity Of 3 
Least Noisy Leads 



603^ 



1 



Set Lead_ctr = 1 




605 



Receive Data For ECG Waveform 



I 



610 



Low-Pass Filter The Data 

1 



-608 



Determine The Dominant Template 
Amongst The Beats To Be Used To 
Construct The Representative Beat 



620 



I 



Create And Store Alignment Template 



Set Beat_ctr = 1 



L630 



Get Next Typical Beat And Normalize 
A Portion Surrounding The QRS 



® 




640 



Fig. 9 A <t> 



23 



EP 0 776 630 A1 




650 





Slide This Data Over The Alignment 
Template While Computing The Minimum 
Area Difference (MAD); Save The Position 
And MAD Value If Minimum For This Beat 




670 


[ Yes 


Increment Lead_ctr 








Fig. 9B 



qso 



Correct Each Beat Time With The Median Value Of 
The 3 Alignment Times Just Determined And Store 



7T~ 

C Return) ^ 



699 



24 



EP 0 776 630 A1 



Create Representative Beat 
From ECG Waveform 



I 



L700 



Determine Which Beats To Use 



705 



I 



701 



Set Lead Num = 1 




^799 
C Return) 



Display 
Representative 
Heart Beats 



780 



Read Aligned Beats For This 
ECG Waveform 



I 



725 



Set Time_ptr = 0 



720 



Get 1 Slice Of Magnitudes From 
Beats At Same Relative Time 



I 




Sort Magnitude s At This Time Slice 
^730 



Fig. 10A 




25 



EP 0 776 630 A1 





7:35 




Discard Predetermined % Of 
Smallest And Largest Magnitudes 



Average Remaining Magnitudes 



I 



Store In Representative Beat Array 
At Correct Relative Time 



755v 



I 



Increment Time_ptr 



765 




No 



740 



750 



Forward/Reverse LP Filter 
Representative Beat 



775 



Increment 
Lead Num 



t 



1 



store 
Representative 
Beat 



I 



770 



Fig. 10B 



26 



EP 0 776 630 A1 



820 




k800 



Read Representative Beats 



I 



^-801 



Measure QRS Onset 
And QRS Offset 



803 



Set Lead_Num = 1 



5 



805 



810 



Get Representative Beat For 
Lead Lead Num 



815 



i 



Determine Isoelectric Level 



I 



Determine R-Wave Amplitude 



I 



Determine ST Level 



I 



825 



Determine ST Slope 



830 




Fig. 11A 



27 



EP 0 776 630 A1 




835 



Update Measurement 
Confidence Flags 



.850 




855 



880 



Display Measurements 



Fig. 11B 



28 



EP 0 776 630 A1 



( Start ) 



150 
1- 



Receive Representative Beats And 
Measurements From Computing Unit 



Print Representative Beats 
And/Or Measurements 



190 



C End 



99 



Fig. 12 



29 



EP 0 776 630 A1 




CO 



CD 

u. 



> 



30 



EP 0 776 630 A1 




31 



EP 0 776 630 A1 



(O 



D) 

9mm 

u. 




7 



0) 
Q. 

E 

<D 



33 



EP 0 776 630 A1 



Beat 1 



Beat 2 




: : Fig. 17 

ft 



Beat N-1 



I 

Beat N 




34 



EP 0 776 630 A1 




• • • 



EP 0 776 630 A1 



« • • • 



• m t • 



• 0 • 
■ « 



• • • • 

• ■ • 



« • • * 



• • • 

• « • • 



t • 



• I 

• « • 



■ 

* 



• • « 



• • ♦ 
« • • 



■ *» 
• • • 



I 



ST licvcl: -2.5 mm 
ST Slope: 4S mm/f ce 
n[ ST Integral: 0.3 cm2 
> Heart Rate: 112 bpm 




• . a 



»• • • 



a 
• 

■.•as 

• • • 

• • • 

• •* 

• •. 

i a A 

• a . 



■ * • . 



p • « a 



a • 
• •a* 

• a a 



• •• 

• •« 

« • 

• *• 

• •• 

• • 



I . 



Fig. 19 



36 



EP 0 776 630 A1 




European p.ten. EURQPEAN SEARCH REPORT 

OfTtot EP 96 36 8432 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Calccory 



CitMioa af docuawal wHh indicatioii, when appropriate, 
oirdeyml passa(<s 



Relevant 

to daim 



a^ssincATioN of thr 

APPUCATION (IntCU) 



X 
A 



US 4 589 420 A (ADAMS ISABELLE H ET AL) 
20 Hay 1986 

* column 3, line 4 - column 11, line 64; 
tables 1-12 * 

US 4 930 075 A (KORTAS RICARDO G] 29 Hay 
1990 

* abstract * 

* column 1, line 43 - column 3, line 11 * 

* column 3, line 37 - column 8, line 20; 
tables 1-6 * 



1 
2 

2-10 



A61B5/0452 
606F19/e0 



The present seardi report has bcco drawn up for all dains 



TECHNICAL HELDS 
SEARCHED (lBt.a.4) 



A61B 

G06F 



BERLIN 



Di<t of coaiplcilo» of tkm i*mth 

14 February 1997 



Weihs» J 



CATEGORY OF QTED DOCUMENTS 

X : pvticulaily relevant if taken iloM 

Y : panicuUriy relevant if combined with aaotlicr 

tecumcnt of tiM sanM calegpiy 
A : lechnolo^cal backgnMnd 
O : non-written disdositre 
P : intcrmeiliate docuincnt 



T : theory or principle undcriying the lavotion 
E : earlier patent docnmeot, but pnblishcd oa, or 

after the fili&g 4ate 
1) : dDcnment dted in the application 
L : document at«4 for other rcasoni 

dfa : nenber of the same patent faaiily, corrcspoaiiag 
docuiaent 



37 



